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ABSTRACT 


The icing of structures due to prevailing climatic 
conditions can present serious problems. Various methods 
to remove the ice build-up have been suggested. One of 
these is thermal radiation. However, before it can be 
used, the factors governing thermal radiation penetration 
in ice must be considered. A computer model was developed 
to simulate the transfer of collimated thermal radiation 
_ through an absorbing and scattering ice sheet. The model 
utilized a modified two flux representation of thermal 
radiation transfer first developed by Schuster (1905). 

The validity of the model was checked by completing 
two sets of experiments. One set of experiments measured 
amounts of incident radiation absorbed and transmitted 
through various ice samples. In order to then compare 
these experimental results to the model, a second set of 
experiments were completed to measure the amount and distri- 
bution of scattered radiation produced in a thin ice sheet. 
Sufficient evidence was obtained by these experiments to 
verify the model. 

The model was then used to calculate the amounts of 
incident energy that are absorbed, transmitted and reflected 
by an ice sheet for the various scattering characteristics 
of ice. This information provided an indication that thermal 


radiation might be useful as a deicing method. 
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CHAPRTERSE 


INTRODUCTION 


The formation of ice on structures whether on the 
land, on the oceans, or in the air, can present very serious 
problems. When icing has occurred, if no viable method 
was available to prevent or failing that remove the ice 
build up, the result has been the collapse of power lines, 
the loss of ships and the downing of aircraft. Realizing 
the problem is only part of the answer. Next it will be 
necessary to propose and evaluate several methods which can 
be used to alleviate the icing problem. 

Depending upon where the icing has occurred, various 
methods have been suggested and some have been tested under 
field conditions. On the land, the problem of icing has 
been solved by first evaluating conditions under which icing 
occurs. The structure is then designed to withstand the 
predicted ice load or is moved to a location where icing 
would not occur. On the oceans, the problem of ship icing 
has been tackled by again evaluating conditions under which 
icing occurs, however, the problem of preventing or remov- 
THgmt hen TCembUl! d Upeiissmore dif hicultstowsolver ss Various 
researchers (1) have attempted to use foam mats and pneu- 
matic tubes attached to ship surfaces to facilitate ice 


removal. In the air, the aircraft industry (3) has used 


ae, ei 
ee Wests 


; a 7 The 

at no re ioe eave aie Z 

aeninse Urey THs a4 nea, sat fi te. ees 7 

bouran niSaiv on rt. prenes 2na “pihk vod aaotdaee | ; 

| aot hid oh al seni git That 9 tnsyetq oF afdettava ae 
~aeatt “eval te sequttae ait need ee ttveey oto «ath ot ted = Ee 

parsttenk MAYO TES to PArWoh sad bon valde to geet ott ‘ 


od itu 77 Jeoh -vewads aft ty tsoq yew et matdovg att 7 
dais totidh svortton reyaves dteuleve bis o7roqerd of yisezonem: <r - | 
maior oniot nde opetvelts of bee “7 

zderiev ,batws56 gen ghist ald en hitiv jody: qe rbasqed : 
“ely betzey \Aes6- 4Veq oatsi> hers bay eHpoya, peed ‘svine abv) om 
ona @ntch ‘ta hardére ot Singh <4F7n0 ~. 2nerti laos ntoiy 
pets? tot4e-vebou enobi yt yon Pere texby ye iévies oh 
—gild -boes eaFhw od ree ‘weata “zt Seubsunta sat aise 
eatot snot wottero! #2 boven 2t ote heal oor sgaotbiong | 
Rntst ghia to mal dg. sis. s BAD 9. did WO .Awoee don bi 


sie ene, Atala pee Shee Nerina ansd & 


various methods including the use of engine exhaust, 
pneumatic tubes and specially treated surfaces to which 

ice has a low adhesion. The various methods of ice removal 
can be divided into two main categories - thermal and 
mechanical. 

Under these two main categories are a variety of 
methods several of which have been illustrated in the pre- 
ceeding paragraphs. The objective of this thesis will be 
to examine thermal radiation as a thermal method to aid in 
deicing surfaces. When using thermal radiation on an ice 
Sheet, the objective is to allow penetration of the ice by 
the radiation resulting in its absorption at the ice- 
surface interface. By absorbing the heat energy at the 
bonding surface it is hoped that the bond could be broken 
allowing mechanical means of ice removal to be more effec- 
tive. The model developed in this thesis, although it is 
intended primarily for study of ice removal techniques, 
could be used to predict such phenomena as transmission of 
solar radiation through lake ice. 

For any problem involving thermal! radiation penetra- 
tion in an ice sheet, it is necessary to Know the effects 
of the type of thermal radiation to be used and the various 
factors which govern its degree of penetration. Previously 
the approach to this problem has been to gather experi- 
mental data, on certain types of ice, in order to determine 
total extinction coefficients for a specific spectral distri- 


bution of thermal radiation. For example, Thomas (4) deter- 
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mined the effects of density and porosity on the extinction 
coetgacients Forssnow) and@ice. | iihekextinctiohacoefficients 


were used in the Bouger-Lambert Law (i.e., I = I e7 KX) to 


0 
evaluate light transmission. In these studies the roles 
of scattering and absorption were not considered separately 
merely their combined effect. A more generalized approach 
to thermal radiation penetration in ice is required if 
thermal radiation is to be evaluated as a deicing method. 

The generalized approach used will consider the ice 
to behave like a diffusing material. That is, thermal 
radiation entering such a diffusing material is absorbed, 
transmitted and scattered within it. This approach was 
first developed by Schuster (5) to describe the absorption 
and scattering of radiation in what was termed a "foggy" 
atmosphere. The equationswere further modified by 
Silberstein to take into account the passage of collimated 
thermateradiation througn a diffuse material. S.Q. Duntley 
(7) discussed further some other alterations of the 
Schuster model to take into account such factors as 
dependency of scattering and absorption on wavelength, and 
the dependency of scattering and absorption on depth, and 
amounts of diffuse radiation. 

It is now hoped that this model first developed by 
Schuster can now be applied with suitable modifications 
to solve thermal radiation penetration in an ice sheet. 
By comparing the theoretical results with those obtained 


from experimental tests, the model can be evaluated. 
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Information obtained from the model can then be used to 
indicate the possibility of using thermal radiation as 


a deicing method. 


, : ‘a: ed vay Adt™ — 
ie) ; 7 ) a, : 
a) oe 
ot bexu ot pate as Talon At’ m0? bambedits iM pine 
q 7A re 
a6 notice’ tan fri 7 nav, bo Ur k PPD a2 08 Bat oiasts at 
oe ; | = 


vos ow gator 


_ 


CHAPHERRI! 


ANALYTICAL ANALYSIS OF THERMAL RADIATION 
TRANSMISSION IN A DIFFUSE MATERIAL 


2. introduction 

For a collimated beam of thermal radiation of known 
energy striking the surface of an ice sheet, the problem 
is to determine how much of the incident energy is absorbed, 
transmitted and reflected by the ice sheet. To answer this 
a modified analysis of the equations presented by Dunkle 
and Bevens (8) will be applied. These differential equa- 
tions along with the specific boundary conditions applic- 
able to the problem will be used to obtain a set of 
expressions which give the magnitudes of upward and down- 


ward fluxes within the ice sheet as a function of depth. 


2.2 Differential Equations 


Consider an infinitesimal layer of ice located within 
an ice sheet subjected to a beam of collimated thermal 
rodlauiondstnikingtnormalbtosthesuppervstrface.sthoraa 
given wavelength and a given interval of depth from x - dx/2 
to x + dx/2 there exist three fluxes. These fluxes are a 
collimated flux and two diffuse fluxes - one in the upward 
and one in the downward direction caused by the scattering 


of the incident beam. To obtain the differential equations 
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some assumptions about the character of the absorbed and 
scattered radiation within a given layer must be made. 
First, the nature of the diffuse fluxes are such that upon 
leaving the infinitesimal layer they are isotropically 
distributed over a cone of given angle (9). Second, the 
proportion of the collimated beam scattered forward (n) 

is different from that scattered backward (1-n). And 
finally, the absorption and scattering coefficients for 

the collimated and diffuse fluxes are related by a constant 
determined by the cone angle of the diffuse radiation. That 
is, it is possible to relate the various absorption and 


scattering coefficients by 


and 


where K and Queetee the absorption coefficients for the 
collimated and diffuse fluxes respectively and similarly 

We and r are the scattering coefficients. Cy is a para- 
meter varying from 1.0 to 2.0 depending upon the cone angle. 
From the above assumptions three differential equations 

can be formulated (Figure 1) 
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dZ = K Zdx + rZdx = rYdx - (1-n)r Idx (2) 


di = -r Idx - KI, dx (3) 


Bate x (1-n)rglydx KZ dx 


kY dx nfel, dix kIydx Rade 


Figure 1. Radiation Fluxes Within an Ice Layer 
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defined as (re + K). Solving the remaining equations 
for the diffuse fluxes is accomplished by defining u = 
Y + Zand v = Y - Z and then adding and subtracting equa- 


EiOnsetyeang (2)etosoptain 
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dxee -K u + het, (5) 


where oh is defined as (2(r}) + Ks Now by differentiat- 
ing equation (1) by — plus substituting the equation for 


I, the following equation is derived: 
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ae 4 by a - o(2n-1)r as (6) 
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By utilizing equation (5) the differential equation becomes 
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The solution of this equation is of the form 
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where re is defined as ak and Cy plus Cy are constants 


of integration determined by the boundary conditions. From 
ThewexpressiOn TOM Ut 1S possible tovobtain one for Vi by 


using equation (4) to get 


De -(-F)e,e>* + ( 
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By replacing some terms by simpler ones and then add- 
ing and subtracting equations (7) and (8), equations are 
obtained giving the upward and downward fluxes. These 


equations are 
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2.4 Application of Boundary Conditions 


In order to establish the values of the integration 


1] 


constants in equations (9) through (11) the boundary condi- 
tions for the various wavelength intervals must be used. 
There exist three boundary conditions: one due to the 
reflection of the collimated radiation at the upper air-ice 
interface and two because of reflection of the upward and 
downward diffuse fluxes at the upper and lower air-ice 
surfaces respectively. 

The boundary condition resulting from the passage 
of collimated radiation through the upper surface (Figure 
2 yeimplies that Cy = (1-r5)I, where a iiSmcnemcoe tic ent 
of reflection for collimated radiation at the surface and 
Ty is the intensity of incident energy on the ice sheet. 


Therefore equation (11) becomes iv = Clete get 


Figure 2. Collimated Boundary Condition 
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The other two boundary conditions are determined if as 

stated earlier the interaction of the diffuse fluxes with 
the upper and lower surfaces is considered. At the upper 
surface or air-ice interface the only portion of the down- 


ward flux existing is that due to the downward reflection 
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(Figure 3). 
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assumed that the reflection coefficient between the air- 
ice interface and the ice-air interface is approximately 
the same (10), which is true for angles of incidence less 
than 30 degrees, and that Cc; can be used to relate the two 
coefficients. Similarly the reflection of the downward 
flux at the lower surface or ice-air interface plus the 
diffuse reflection of the collimated flux produces the 


final boundary condition (Figure 4). 


Figure 4. Diffuse Flux Boundary Condition 
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For the lower surface the relationship between ers and re 
is given by the expression rp Su Ss eae Als om tees 
assumed that the reflection of the collimated radiation 
is to be totally diffuse. From the above developed boundary 
conditions it is now possible to evaluate the integration 
constants for the diffuse fluxes. 

If the equations derived earlier at a given wavelength 


are evaluated at x = 0 and x = b the diffuse flux equations 
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Applying these flux equations to the boundary conditions 


two equations are obtained 
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Solving the above equations expressions for the two con- 


stants are 
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CHAPTER 1it 


COMPUTER SIMULATION 


3.1 Computer Model 


Using the equations developed in the preceding , 
chapter, the radiation fluxes at a given wavelength and 
dapchieinmdneicessncetacan be found. s) Inesesimuxesmalong 
with the absorption coefficients (Figure 5) - obtained 
from data presented by Goodrich (11) - are used to calculate 
the portion of incident energy that is absorbed by an ice 
Sheet. The incident energy transmitted and reflected are 
evaluated using the radiation flux equations at the sur- 
faces of an tce sheet. 

To numerically evaluate the iene aaa absorbed, 
the ice sheet is divided into a series of finite layers. 
The radiation fluxes existing in each layer are determined 
at the center of the layer and assumed constant over the 
remainder. Using the absorption coefficients, the energy 
absorbed within each ice layer is then found by multiplying 
the fluxes by these coefficients. Since this absorbed 
energy is evaluated at a given wavelength, the total radia- 
tion absorbed within each ice layer must be numerically 
integrated over the total spectrum. This is done by the 
trapezoidal rule. The total energy absorbed is then deter- 


mined by adding the energy absorbed within each layer. 
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The same trapezoidal rule can be used to find the 
total incident energy transmitted and reflected. However, 
first the amounts of the transmitted and reflected radia- 
tion must be determined as a function of wavelength. To 
do this, the boundary conditions and the flux equations at 
the upper and lower surfaces must be used. The boundary 
condition are Gay, + Obata it. for the transmitted radia- 
tion and (lena) 2 + rol for the reflected radiation. To 
compare the effects produced by the different incident 
energy spectra used, the absorbed, transmitted and reflected 
energies are normalized with respect to the total incident 
energy. The total incident energy is obtained by numeric- 
ally integrating over the total energy spectrum using the 
same method employed to determine the transmitted and 
reflected energies. By summing up the normalized values 
of the absorbed, transmitted and reflected energies, a 
check on the numerical accuracy of the computer model is 
obtained. The summation indicated that for a 20 cm thick 
ice sheet the maximum error was two percent. This particu- 
lar ice sheet had been divided into 50 finite layers and 
had a scattering coefficient of 1.0/cm. 

Due to the nature of the flux equations used, however, 
several special cases arise in computing the amounts of 


energy absorbed, transmitted and reflected. The first 


problem that occurs is when o is equal to @ and the term 


a appearing in the diffuse flux equations cannot be 
(a) - 
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evaluated. By setting a = = and then using the relation- 


Ship between ne and re plus K and a the critical condi- 
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Therefore, for a particular degree of scattering there 
exists a critical absorption coefficient for which the 
diffuse fluxes are impossible to evaluate. This numerical 
problem results from the form of the particular solution 
chosen to solve the diffuse flux equations. To overcome 

Crem noblenptnesvalvesOnekals resetasiightivaewnenmitaf alls 
on the critical value. Since this small change is insigni- 
ficant considering the uncertainty in the absorption co- 

Gt ficients, (11) no important error is» produced. “The next 
special case occurs when the absorption coefficients are 50 
times greater than those for scattering at a given wavelength. 
When this situation occurs, it is assumed that the passage 
of the collimated radiation obeys the Bouger-Lambert law 

and that only absorption need be considered. The absorption 
is then determined by evaluating the difference between the 
collimated flux at each surface of the finite ice layers. 
Finally, numerical problems occur in the model when an ice 
SNECCtmISMODETCallYaVervaetnick ya chatmls, one and Ch is very 
large. When this occurs the diffuse flux integration 
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Utilizing the general case and the various special 
cases for radiation penetration in an ice sheet a series 


of computer outputs are obtained. 
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The major parameters which affect the results are the 


temperature of the blackbody source, the degree of non- 


isotropic scattering (n), the scattering coefficients and 
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the diffuse cone angle constant (cy). Two series of curves 
are generated. One series Appendix B gives the behavior 
of the absorbed, transmitted, and reflected energies for 
various combinations of the major parameters. The second 
series Appendix C gives similar sets of parametric curves 
for the rates of energy generation in the ice sheet due to 
absorption of the radiation fluxes. From these parametric 
curves the effect of varying the major parameters are 
evaluated. 

A decrease in blackbody source temperature caused the 
radiant energy to shift to the longer wavelengths. Because 
of the higher absorption coefficients that exist at the 


longer wavelengths (Figure 5) . this produces an increased 


absorption, decreased Benirsiieiaion and decreased reflection 

of the incident energy. As an example for a 12 centimeter 

thick ice sheet with no scattering, for a 3000°K source 

temperature 38 percent of the energy is transmitted while 

for the 5780°K source temperature 72 percent is transmitted. 
The effect of scattering is quite important as far 

as the transmitted and reflected energies are concerned. 

In general, the effects of scattering on absorption are 

minimal. Increased scattering produces a decrease in trans- 

mission and an increase in reflection. This effect is 

particularly important for the high temperature source. 

Due ‘to the ‘higher absorption coefficients ithat sexist at ithe 

wavelengths where the lower temperature source is the 


strongest, the effects of scattering on the transmitted and 
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reflected fluxes are less important. As would be expected 
these effects of scattering are not linear. That is, 
because of the exponentials in the flux equations, equal 
increments of scattering coefficients do not cause equal 
changes in the transmitted and reflected energies. The 
effect of scattering on absorption is minimal as stated 
earlier. However, changes in degree of forward scattering 
and scattering itself do produce some interesting phenomena. 
For example (5780°K and n = 0.90) an increase in scattering 
causes an increase in absorption for ice sheet thicknesses 
less than 12 centimeters. While for thicknesses greater 
than 12 centimeters, the absorption is decreased with 
increasing scattering. This behavior can be traced to two 
opposing effects caused by scattering. First, an increase 
in scattering produces a greater amount of diffuse radiation 
in the ice sheet. As a result of this, more incident energy 
is absorbed since diffuse radiation has effectively higher 
absorption coefficients compared to collimated radiation. 
Second, as more incident energy is scattered, a higher por- 
tion of energy is reflected out of the ice sheet. These 
two effects then counteract each other causing the observed 
behavior. 

The degree of forward scatter as opposed to backward 
Scaveerealsouhas ani eLrect.enm [US moOStusioniiicant.errege 
is that which it has on the incident energy absorbed. By 
increasing the degree of forward scatter (Appendix B, pages 


93 to 98) the crossover of the absorption curves for the 
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various scattering coefficients is delayed. This delay is 
a result of less incident energy being reflected back out 
of the ice sheet. The other effects of increasing forward 
scatter are an increase in transmission. 

To determine the effect of the diffuse cone angle 
constant (cy) a comparison is made between the generation 
curves (Appendix C, pages 108 and 109). These are produced 
by varying the value of Cy from one to two at a given 
degree of scattering coefficients. This comparison indic- 
ates increasing Cy has a similar effect to a slight increase 
in scattering. Now, a further examination of thermal radia- 
tion penetration in an ice sheet will require input from 
experimental work. By completing these experiments repre- 


sentative values of the major parameters can be chosen. 
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CHAPTER IV 


EXPERIMENTAL PROCEDURE AND RESULTS 


4.1 Introduction 

The two diffuse flux model with collimated flux super- 
imposed will be tested for its validity by using two differ- 
ent types of experiments. The first type of experiment 
will establish the accuracy of the predictions for the 
amounts of incident radiation absorbed and transmitted; 
two different sources (Plate 1 and Plate 2) providing the 
necessary incident radiation. The second type of experiment 
will then be used to test assumptions about the character 
of the scattered radiation by determining what proportion 
is scattered forward and backward along with the scattered 
energy distributions. Finally, by comparing the informa- 
tion obtained from the experiments with that obtained from 
the computer model, some conclusions can be drawn about 


radiation transfer in ice. 


4.2 Experimental Apparatus 


4.2.1 Temperature Controlled Chamber and Ice Calorimeter 


To determine the amounts of incident energy absorbed 
and transmitted, an ice sample is placed in a temperature 


controlled chamber (Figure 6 and Plate 2). This chamber is 
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maintained in the temperature range from -8°C to 0°C for 
the various absorption and transmission measurements. The 
ice sample has embedded in it three thermocouples which 
are frozen into the sample. Located within the chamber 

is an ice calorimeter (Figure 7 and Plate 3) into which 
the ice sample is placed. Surrounding the sides of the 
ice sample is a metal reflecting surface. Its purpose is 
to simulate, as far as the incident energy is concerned, a 
one-dimensional ice sheet. To allow the incident energy 
to strike the ice sample and then measure the portion of 
energy transmitted a flap is located on either end of the 
calorimeter. Two openings are also provided on either end 
of the temperature controlled chamber. Located on the 
front of the chamber is a two millimeter water filter 
(Figure 8). Its purpose is to eliminate the longer wave- 
lengths contained within. the incident energy spectra. The 
other end of the chamber is covered by a filap: By removing 
the flap, a thermal radiometer can be periodically exposed 


to the transmitted radiation. 


An2. 2  Distrubutadon, 0feccatteredukRadg ation 

The second experimental apparatus was employed to 
measure the distribution of scattered radiation from an 
ice sheet of approximately one centimeter in thickness. 
Basically it consists of a monocromatic beam of light 
(0.63 u) produced with a laser and directed toward the ice 


sheet through a series of lenses. A more complete descrip- 
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tion of the laser used as well as other instrumentation 
employed in the experiments is presented in Appendix D. 
Supporting the ice is a mount to which a quantum radio- 
meter has been attached. The radiometer is attached in 
such a way that it is able to rotate through 60 degrees on 
e1thereside of theslaser beam (Plate 4) Using the appana- 
tus, data was collected from which the angular distribution 


of scattered light was subsequently determined. 


Ave? om Seby Une tic ata Olee = Thermocouples and Radiometers 


The copper constantan thermocouples placed in the 
ice sample are surrounded by a stainless steel shield 
(0.083 centimeters in outside diameter) and calibrated 
to within 0.05°C using a platinum resistance thermometer 
traceable to the National Bureau of Standards. These thermo- 
couples are held in a plexiglass mount which is frozen into 
the chosen ice sample. This procedure allows the thermo- 
couples to be embedded in the ice at 0.8 centimeter (thermo- 
couple #1), 6.5 centimeters (#2) and 11.2 centimeters (#3) 
inerehbauionetouthe Upperesunrace. = Pigures9 islustrates 
how the thermocouples are wired together to obtain the 
necessary temperature reading. These readings were recorded 
on a Hewlett Packard strip chart recorder. 

As stated earlier two radiometers were employed to 
measure radiation intensities. To measure the intensity 
of the incident beam and the transmitted energy a Thermal 


Radiometer was used. This particular radiometer had a flat 
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response for a range of wavelengths from 0.25 microns to 
15.0 microns. To calibrate the radiometer a blackbody 
source was utilized (Figure 10 to Figure 11). For the 
measurement of the scattered radiation distributions a 
Quantum Radiometer was used. Although this radiometer does 
not have a flat response similar to the other probe, only 
One wavelength is involved in the measurements, therefore, 


no corrections are required to the probe's output. 


4.2.4 Ice Mould 

To produce uniformly sized ice samples, a plexiglass 
mould 15 centimeters square and 12 centimeters deep was 
constructed. On the bottom provision was made to attach 
a copper plate cooled by a chilled water-alcohol mixture 
(Figure 12). Depending on the desired ice sample, a plexi- 


glass cover could also be attached on the top of the mould. 


4.3 Experimental Procedure 


ATS We eADSOLDLLON. A ansSMisSdonsand Reflection Measurements 


The experimental measurements of the absorbed, trans- 
mitted and reflected portions of the collimated incident 
energy were obtained for three ice samples and two radiation 
sources. The two sources used were a 650 watt tungsten- 
halogen projector lamp and a Metro-Lite solar simulator. 
Radiation intensities for the two sources with and without 


the water filter are compared to those of a blackbody 
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(Figure 13 to Figure 16). These intensities and source 
temperatures having been obtained from Gilpin (16) are 
found to be comparable to those from an equivalent black- 
body source. Before a given ice sample is exposed to the 
radiation, it is allowed to reach an equilibrium tempera- 
ture with its surroundings. The sample is then exposed to 
the radiation source. If the tungsten lamp is used the 
time of exposure is 10 minutes and if the solar simulator 
is used then exposure is for 30 minutes. The flaps on the 
calorimeter are then lowered and the gradient of tempera- 
ture within the ice is allowed to disappear. In the case 
of the tungsten lamp the time required is approximately 60 
minutes. Similarly for the solar simulator the time 
required is approximately 40 minutes. One problem en- 
countered in the experiments is that heat is transferred 
from the ice sample to the chamber during the experiments. 
To determine the amount of this heat transfer several tests 
were conducted in which temperature changes within the ice 
were measured as a function of time when no incident radia- 
tion was applied. The tests were conducted with the flaps 
open and also with the flaps closed. 

Using the heat loss rates (Figure 17 and Figure 18) 
plus the thermocouple error (0.05°C), estimates of total 
error were made for the absorption measurement. The error 
produced by the thermocouples results in an uncertainty in 
the absorption measurements. However, that error produced 


by heat transfer between the ice sample and the chamber can 
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be approximately corrected for. Utilizing the average 
chamber and ice temperatures measured during each test, 
the heat transfer is estimated from Figures 17 and 18. 
These heat transfer rates are calculated separately for 
the flaps open and the flaps closed. The results are then 
used to estimate the apparent increase or decrease in the 
absorption measurements. The reliability of the absorption 
measurements was then ascertained by making three runs of 
any given ice sample and radiation source (Figure 19 and 
Figure 20). Each curve presented being the mean value of 
three individual runs. 

During each run with the various ice samples, read- 
ings from each of the thermocouples are taken every 5 
minutes when the flaps are open and every 10 minutes when 
the flaps are closed. The radiometer was used at the same 
time to measure the transmitted radiation. However, except 
for a clear ice sample (no scatter), this particular mea- 
surement was only taken after the final run of any given 
ice sample. This procedure was required since the thermal 
radiometer did not function when placed within the tempera- 
ture controlled chamber if it was running. The reason for 
this change was that when the transmitted radiation was 
totally diffuse the measurements had to be taken as close 
‘to the ice sample as possible. Since the chamber was shut 
down for the transmission measurements, it was possible to 
take several measurements of the reflected energy before 


one of the ice samples deteriorated. By converting the 
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raw data, the values of the absorbed, transmitted and 
reflected energies can be obtained in percent of the inci- 


dent energy. 


An3e-72 scattering Distributions 


By mounting a thin ice sheet in the laser apparatus, 
measurements of the energy distributions occurring due to 
scattering were taken with the quantum radiometer. The 
ice sheets were cut from the ice samples used in the first 
set of experiments. Each ice sheet had. both its surfaces 
made parallel by placing the sheet on a heated copper plate. 
Total energies striking the radiometer's detector surface 
are recorded for angles measured from the direction of 
propagation. The angles are in 10 degree intervals from 
60 degrees to 300 degrees and from 120 degrees to 240 
degrees. Using the solid angle subtended by the detector, 


the intensities of scattered radiation can be calculated. 


io Ome COmoalipi @'S 


In order to obtain the required ice samples with 
varying scattering characteristics, the dissolved gases 
present in water were entrapped during ice formation. By 
varying the rate of freezing or the initial gas content 
of the water, a desired ice sample is produced. The above 
method is used to produce three ice samples. Of the samples, 
one was clear and two contained gas pockets acting as 


scattering centers. In one scattering sample the gas cavi- 
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ties were a order of magnitude larger than the other. This 
particular sample was produced by rapid cooling of water 
containing dissolved carbon dioxide giving the ice sample 

a white appearance (Tests #7 to #12). As for the remaining 
scattering sample, it was produced by rapid freezing of 
heated tap water obtained from the City's water supply. 
This ice sample appeared to be gray in color similar to 


Dive CEMULCS CSE smLOMTIG ). 


4.4 Conversion of Raw Data 


4.4.1 Calculation of Absorbed Energy 


Let T be the measured rise in ice temperature 


ri sie 
from the beginning to end of each run. Therefore, the 
energy absorbed by the ice is ErceE = Mico lise where my is 
equal to the mass of ice and Cy is the specific heat. Some 
of the collimated energy that is absorbed by the ice is 
conducted into a stainless steel support surrounding the 
Tecesanp lees SASSUMING etic SUDO tal Sedt athemsamem Lempe na ture 
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rise 
io. cal where Ch of thesstainless steel ts 02 lecal/gqneG and 
the mass of metal used is 554 grams. The percent of energy 
absorbed can now be determined if the total collimated 
energy. striking,the ice surface is found. .Jo,do this, sthe 


ice surface receiving the energy was assumed to be 225 


square centimeters. Therefore, the energy incident on the 
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measured intensity of radiation at the ice surface obtained 
by using the thermal radiometer. The fraction of incident 
+ (0.06)10° T 


m.c é 
then absorbed is approximately { Ip Anes | 


(3.23)10° es 


an 422 elransmitted Energy 


Two types of transmitted radiation are measured. The 
clear ice sample, since there was no scattering involved, 
produced a collimated beam while the two scattering samples 
produced a totally diffuse beam. When the radiation strik- 
ing the thermal radiometer is diffuse, only a portion of 
the radiation is detected by the radiometer. This propor- 
tion is found by considering the internal geometry of the 


radiometer's detector. For this particular thermal radio- 


meter employed, the portion of energy detected is 1.94 x 10° 


The fraction of collimated energy transmitted is found by 
comparing the radiometer measurements, corrected if the 
radiation is diffuse, with those taken for the intensity 


of the unattenuated collimated beam. 
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The quantum radiometer measurements of scattered 
radiation indicate the amounts of energy intercepted by 
the detector. To convert these measurements to radiation 
intensities, the solid angle subtended by the radiometer 
must be found. Using the detector area plus the distance 


between the centre of the thin ice sheet and the detector, 


Z 


an 7 a a“ a 
g 0. ee ¢ 


a? 


- 


mnt jae Aone 


- ae Bs ary 
‘ ; 7 a - 
| iy +6 ‘eit _— 
bas tapos 42 12 Bot She of iota tos" YO VO eatery. Be 
: he —) 
pnehtonh to wotsonrt. edt _ Sopenet be rein 
~ ‘ . 
raed aOi0) +o ee . 7 
{ tot he = . Oy Vie BS “908 ba fi 
1) OU PS,E) ae ¢ 
} - 
pe ees ; 
{ ey 
vovena .h MAD Coe i. 
AG, me oe ae traced f 
i27 “6 notdotban bage finan eid '% basin citer ra 
OVsTi rey ogy a- OF 244 ie ty ae ‘¢? a Kg? ¢ aut spats 
; - na ‘ iis oe 
. s4i5ge O i+ ah paw hateat (loo s — Me 
VT r3 ip ‘ J | fm iD x i ai 07 & re saat i. 
: es baa dl 
oF ; 0) be fi OND) Dea sata wit 
d ame then. 3 vel a5 ol oF 
P a? é aa of i 
r AA SA bebe aD, CH Ago ee 
; : an } i) 
rusian 2 109755 ob 
979 R9D VTS 4 YT. eT See Syl , 4 bay, 
_ ’ ‘ 4 4 Om. a4 Dé 
5 bada ine yuyernens bed eae Too? 
r * ia > 
; ho i. eden $ot smoteay 
, , i* eh eye 
a | ‘s 
v oth { fax so ‘920 ig ha ees 
Vi. 
at bf che A fi [gas 
i 
y (Soy ea 8 0 a2Homarvagee 30 io FOS git sel 
- 7 
ea beaiaat Eas eis i Se cee 
: 7 a t 
te pie 
ne f mW sve. ies 
Jeera th pe saad 1 srenwinge 
ee 


45] 


the solid angle was found to be 0.0676 steradians. The 
scattered radiation intensities are found by dividing the 
radiometer reading by the detector area and the solid angle. 
These radiation intensities are then normalized with respect 
to the incident energy striking the ice. 

To determine the fraction of the collimated beam that 
is scattered, comparison between the energy incident on the 
detector when it was set at zero degrees with no ice sheet 
present and the output of the radiometer for the attenu- 
ated laser beam is made. The Bouger-Lambert law is then 
used to determine the scattering coefficient. By assuming 


absorption is small compared to scattering (true for 0.63 u 


-Y_X 
in ice), the law becomes I = Ie oS Solving this equation 


gives the value of the scattering coefficient ee 
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The data obtained from the absorption and transmission 
experiments are presented in Tables I and II while those 
for the distribution of scattered radiation are presented 
imeeigure.2| sands bable wliie SErom thesexperiments isevera| 
observations can be drawn about radiation penetration in a 
12 centimeter thick ice sample. The first of these involves 
the effects of scattering on absorption, transmission and 
reflection of the collimated incident radiation. Comparing 
thea data obtained for clear ice) to thatewrorm the remaining 
ice samples, shows an increase absorption, a decrease in 


transmission and for the limited data available an increase 
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EXPERIMENTAL RESULTS 
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NORMALIZED SCATTERING DISTRIBUTIONS 
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Ice Sample 


Weer ae 2 «8a el Yea e9.0/.cm 
(Degrees) (1/Steredian) (1/Steredian) 
0 7750 20 a0 
10 0.475 0.629 
20 OF 3 OrrZ 19 
30 ORS | 0.149 
40 0.100 0.106 
50 0.082 0.079 
60 0.060 0.061 
120 0.019 0; 201313 
130 0.020 0.036 
140 0.020 0.036 
150 OF 0:20 0.036 
160 = 0.036 
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in reflection. However, changes in absorption with 
increased scattering are less than those for the trans- 
mitted energy. 

The estimates of experimental error for absorption 
measurements indicate that heat transfer from the chamber 
to the ice sample is significant. Since the average chamber 
temperatures were observed to be higher than the average 
ice temperatures, the heat transfer must cause an apparent 
increase in the absorption measurements. A more realistic 
value of the actual energy absorbed would, therefore, be 
obtained by subtracting the increase in absorption due to 
heat transfer from the absorption measurements based on 
the total temperature rise. 

The final group of observations, about radiation pene- 
trating in the ice sheet, are obtained from the scattering 
distributions. In general the radiation is scattered pre- 
dominantly in the forward direction. When comparisons are 
then made between the forward and backward components of 
the scattered radiation, the scattered radiation distribu- 
tion forward is dependent on the angle from the direction 
of propagation while the backward scattered radiation is 
isotropically distributed. The size of the scattering 
centers also has an effect on the scattering distributions 
with the smaller scattering centers producing a greater 
portion of backward scattered radiation. Since for the two 
scattering ice samples, the amount of radiation scattered 


from the collimated beam is the same, the reduction in size 
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of the scattering centers appears to reduce the fraction 
ome energyescatiteredstorward (1.e..9n).0 Usingstnesabove 
observations comparisons can be made with those results 


generated by the computer model. 


4.6 Experimental and Theoretical Results Comparison 


The experimental data must now be compared with the 
theoretical computer model. Since the measurement of the 
angular distribution of scattered energy indicates a high 
forward scatter, the computer model with 90 percent forward 
scatter (n = 0.90) is used for these comparison purposes. 
In general, increased absorption is found with scattering 
(12 centimeter ice sheet) and transmission measurements are 
found to agree with the computer model (Figure 22 and 
Figure 23). Also the limited data available for the 
reflected radiation has good agreement with theory. 

Any error that does exist can be traced to a variety 
of reasons. Some of these are related to the radiation 
sources used. Other sources of error are the possibility 
of non-uniform ice samples or the methods used to measure 
absorption and transmission. Comparisons of the radiation 
sources used for the experiments and the blackbody source 
with water filter are in close agreement (Figure 13 to 
Figure 16). From observations made during the experiments, 
the tungsten lamp is the steadier of the two sources used. 
Except for the clear ice sample, this did not seem to affect 


the transmission measurements which were all within 5 per- 
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cent of the theory. However, when the clear ice Sample was 
exposed to the solar simulator, the measured transmission 
was 15 percent below that predicted by the computer model. 
This larger error is most likely due to the non-uniformities 
in the ice sample. Since these affect the higher tempera- 
ture sources to a greater extent, the error is larger. 

These non-uniformities in the ice samples, however, do not 
affect the other transmission measurements. This is prob- 
ably due to such a large degree of scattering that all radia- 
tion fluxes in the ice are totally diffuse. This point is 
born out by observing the transmitted radiation in the 
scattering ice samples which is totally diffuse. 

Although the transmission measurements reasonably 
match the model, those for absorption do not. The error 
for the tungsten lamp source is 10 to 20 percent and that 
for the solar simulator is 20 to 30 percent based on total 
incident energy. Since the well accepted Bouger-Lambert 
law is used by the computer model for the clear ice cases, 
there must be some inherent error in the experiment itself. 
This error could be produced by two factors. One is the 
difficulty in simulating the one-dimensional ice sheet and 
the other is the heat transfer from the chamber to the ice 
sample. The lack of one-dimensional ice sheet affected 
absorption more than the transmission, since the trans- 
mission measurements are taken along the center line of the 
ice sample. Because of this, the radiation transmitted is 


not affected by the edge effects due to the presence of the 
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reflecting surface. 

As described previously, the effect of heat transfer 
from the chamber to the ice sample was corrected for in the 
calculation of the energy absorbed by the ice. Having done 
this, the difference between the measured absorption and 
the absorption determined by the computer model was reduced. 
The resulting error for the tungsten source was reduced to 
5 to 15 percent and that for the solar simulator to 15 to 
25 percent. This difference between the experimental and 
theoretical absorption of ice, along with the large vari- 
ance in the calculated absorption for the individual experi- 
mental runs, indicates that reliable measurement of absorp- 
tion is extremely difficult. The problem being to expose 
an ice sample to thermal radiation and at the same time 
isolate the sample from its surroundings. However, signi- 
ficant agreement is available from the transmission and 
reflected energy measurements to verify the model chosen. 
Utilizing the information obtained from experiments and the 
computer model, the questions posed in the Introduction can 


be answered. 
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CHAPTER V 


CONCLUSIONS 


To determine the possibility of using thermal radia- 
tion as a deicing method, a knowledge of the behavior of 
thermal radiation in an ice sheet must first be obtained. 

To do this a model of radiative transfer in an ice sheet 

was first developed. The model assumed that ice behaved 
like a diffuse (turbid) material). Experiments were then 
carried out to evaluate validity of the model. From the 
data which was gathered, sufficient evidence was obtained 

to suggest that the model was reasonably accurate. Cer- 
tainly within the accuracy necessary for practical purposes, 
ChismisaoVUeow bnat 1S.. 1neti1eld situations, mene uncertain cy 
about the scattering properties of ice would probably be 
greater than any error in the model itself. 

The computer model is then used to provide the informa- 
tion required to evaluate thermal radiation as a deicing 
method. It generated for various combinations of radiation 
source temperatures and ice sheet scattering characteris- 
tics, the amounts of incident energy that were transmitted. 
in general, it was found that the higher Sources temperatures 
produced the larger amounts of transmitted energy. However, 


asmsource temperatures are increascdms Lnese) fects) 01 scalter= 


ing became more pronounced. 
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To summarize the effect of scattering on transmission, 
consider the mean free path length of a photon before 
scattering of it occurs. The mean free path length is 
expressedias thesreciprocal: Of ethe scabteringmcoetiicient 
(12). For example, if the scattering coefficient is 0.5/cm, 
the mean free path length of the photon is two centimeters. 
It would be expected that, the transmission of thermal 
radiation is not hindered greatly for ice sheet thicknesses 
less than one mean free path length. This is in fact born 
out by the computer model predictions for the 5780°K radia- 
tion source. The reduction in transmission (Appendix B, 
page 85) was approximately 15 percent lower than the clear 
ice transmission. A further examination of scattering and 
its effect on transmission can be obtained by considering 
for what ice sheet thicknesses scattering causes a 50 per- 
cent reduction in transmission. For the same 5780°K source 
temperature (Figure 23, page 63) the ice sheet must be 
about five mean free path lengths ow thick. These con- 
clusions can now be utilized in the consideration of thermal 
radiation as a deicing method. 

In general, the computer results show that thermal 
radiation transmission through thin ice sheets (less than 
2 cm) is significant regardless of scattering when a high 
temperature (5780°K) source is employed. This in turn 
indicates that thermal radiation would only prove valuable 


if the ice build-up was not excessive. Since certain 


objects such as antenna and guy wires cannot withstand a 
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Targe ice build-up, thermal radiation could be a viable 


method for deicing in those circumstances. 
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APPENDIX A 


COMPUTER PROGRAM 
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MAIN PROGRAM---USED TO PLOT DATA OBTAINED BY VARYING 


THE SCATTERING COEFFICIENT FOR THE 
DIFFERENT BANGES OF ICE THICKNESS'S, 


INTEGER TEMP 

DIMENSION ALPA(20) , TRANS (10) ,REFLEC(10) , ATTEN (10) 
7TH {10) , XX (10) ,TSTORE (3,5, 10) 

EQRMADIPANOSE XLS , 2R7 TOP 2 FEA. 2) 

FORMAT (29A4) 
READ (8,101) RS,TEMP,ITH,CI 

IOU=7 

poses 18l=1,5 

RS1=RS* (I-1.9) 

CALL ICE(RS1,TEMP,ITH,CI,TRANS, REFLEC, ATTEN,TH) 

AS=ITH/5.0 

pO 132)3=1710 

TSTORE (1,1,J) =TRANS (J) 

TSTORE (2,1,J) =REFLEC (J) 

TSTORE (3,1,J3) =ATTEN (J) 

CONTINUE 

CONTINUE 


PLOTTING OF DATA 
DOse cen =1 55 

READ (8,201) (ALFA(II) ,11=1,20) 
Pemsier=745 

NF=L 

IF (L.EQ.5.AND.K.EQ.3) NF=-5 
DO 35° M=1,10 
XX (M) =TSTORE (K,L, 4M) 
CONTINUE 

GAELECGPE (£0 ,XSREE RIOT NE, 1 PA eRirOne 
PAS 7535570507027 SS0FALEA PION) 
CONTINUE 
CONTINUE 
STOP 
END 


SUBROUTINE ICE---CALLED BY MAIN PROGRAM TO EVALUATE 


THE TRANSMITTANCE, REFLECTANCE, AND 
ABSORBTANCE FOR TEN DIFFERENT ICE 

SHE OU) CHUCK NESS'Ss ghee USB OULINE 
CONTAINS WITHIN IT FOUR OTHER 
SUBROUTINES USED TO CALCULATE INCIDENT 
RADIATION, ESTIMATE PROGRAM ACCURACY, 
AND HANDLE SPECIAL CASES OCCURING IN 
ICE. 
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SUBROUTINE ICE(PVR,TEMP,ITH,CI, TRANS 


Q ,REFLEC,ATTEN,TH) 


INTEGER %,TEMP 
REAL KAPPA, MEW,MMEW,IN,N 


COMMON WL,IN,T,X (50) ,CABS (59) 
DIMENSION ABS(50,28),TABS (50) , EXC (50) , TOT (28) 


Q ,TRANS (10) ,REFLEC(10) , ATTEN (19) ,TH (10) 
C ICE THICKNESS AND SOURCE TEMP 
Cc 
T=TEMP 
N=0.99 
RC=0,92 
RC1=RC 
DELWL=0.1*1.9E-04 
Honea PALL=1)76 
TH (LL) = (ITH¥LL) /10.0 
DELX=TH (LL) /50 
X (1) =TH (LL) /100 
DOMC29E=0950 
X {L) =X (L-1) +DELX 
02 CONTINUE 
TOTIN=0.0 
ATTEN (LL) =0.9 
€ 
E 
100 FORMAT (F4.2,4X,E7.1) 
200 FORMAT('1',15X,"RADIATION ATTENUATION IN AN ICE 
Q LAYER‘) 
300 FORMAT (? *,10X,* ATTENUATION! ,2X,F6.3) 
319 FORMAT(! ',10X,*RADIATION SUMATION', 2X,F6, 3) 
00 FORMAT ('O',7X,F7.3,17X,F7.5,10X,F8.5) 
506 FORMAT(’? ',5X,'DEPTH',15X,'ABSORTION® 
Q,19X, EXTINCTION COEFFCIENT') 
700 FORMAT ('0!,25X,12) 
800 FORMAT(' *,2X%,5(E15.7,2X)) 
900 FORMAT ('0',20X,E15.7) 
920 FORMAT (*0!,2X,4 (E15.7, 3X) ) 
1000 FORMAT('O',10X,*REFLECTED RADIATION TOP 
Q SURFACE',2X,F6.3,10X,"TRANSMITTED RADIATION", 2X 
O 87603) 
1100 FORMAT('0',10X,'TEMPERATURE=',F8.2,10X,'THICKNESS 
QO =",P6.2) 
Cc ATTENUATION EVALUATION 
1200 FORMAT('0',15X,*RS=',F5. 3) 
c 


DO 03 Z=1,28 
READ(5,100) WL,PVK 


WL=WL*1.0E-04 


74 


as a 
rus pyle ae avays 


i 
1) a ws ‘ 
~~? 


ok 
oun et. tc ‘ 
VE%, (OPE SEND TET ON) 
it, (OR RR RIN 407} 
INST S38 
i +9 
J 
. re (a4 
Chika Zayere tt) x, 
(; 4 a1 oJ 50 rele 
| range (hay E 
ia yng THO 
* § alert BT< 
ne eee 
i” 
: ' 
® x0 Ae. ¥ 
7 
’ . {PF a. ae2 ar 
; al i si ' ATA 4 t yar. a i Nae 
s af r 


ote fe: ‘yj Mer raiiiatee waht 


IM ys, MIRA RG hie Be ft ) 
; oy eit cr T¢ 5 BAD. am J at bt, 
a SERRE EES poke is ees 
(7 LT Sd Tt 
A - 
; 7 if 4 7 A ne tat 


ess 

Li 

zi Ad + salah eeT> 2. 
ss. utes te aaerrienne ghar, 

a ¥ eg ey? ioe 

: i PEAY WS fy eG hee 


a 


23 
28 


15 


CALL ENERGY 
GOsTON28 
PVK=0.95*PVK 
PVK1=PVK*CI 
PVR1= (1. 0-N) *PVR¥CI 
IN=IN*EXP (-PVK*C.20) 
SIGMA=PVR+PVK 
TOT (Z) =IN*DELWL 
IF (PVR. EQ.0.0.OR.PVK.GT. (PYR*50)) GO TO C4 
SIGMA 1= (2*PVR1) +PVK1 
KAPPA=- (SIGMA1+SIGMA* (2. 0*N-1.9)) * (1-RC) *LN* (PVR/2, 0) 
ROW= (SIGMA1*PVK1) ¥*0.5 
IF ( (ROW/SIGMA) .GT. 0.995. AND. (ROW/SIGMA) «LT. 1.005) 
QO ScorTou29 
PI=KAPPA/ ( (SIGMA**2) - (ROW**2)) 
BETA=((2.9*N-1.0) * (PVR/2.0) * (1-RC) ¥IN) 
IF ( (ROW*TH (LL)) .GT.50.0R. (SIGMA*TH (LL) ) 
OQ SeTs50) ecotzoigs 
TAH=1-(ROW/SIGMA1) 
SI=1+ (ROW/SIGM4A1) 
MEW=1+(SIGMA/SIGMA1) 
MMEW=1- (SIGMA/SIGMA1) 
CONSTANT EVALUATION 
R= ((RC*CI) ¥PI*¥MMEW) - (BETA/SIGMA1) *(1.0+ (RC*CTI) ) 
Q -(PI*MEW) 
RR=RC 1* (1-RC) ¥IN+RC1*¥CI*MEW*PI-MMEW* PI 
Q +(BETA/SIGMA1) * (1. 0+ (RC1*CZ)) 
A= (SI-(RC1*CI) *TAW) 
AA=(TAW- (RC1*CTI) *STI) 
B= (SI- (RC*CI) *TAW) 
BB= (TAW- (RC¥CI) *ST) 
DIV=((AXB*¥EXP (ROW*TH (LL) )) — (AA*BB/EXP (ROW*TH (LL) ))) 
CON=- ( (R¥AA/EXP (+ROW*TH (LL) ) )- (RR*¥B/EXP (+S IGMA*TH (LL) 
OF) /DLIV 
CCON=( (R¥A*EXP (ROW*TH (LL) )) - (RR#¥BB/EXP (+SIGMA*TH (LL) 


Q )))/7DIV 
DETERMINE ATTENUATION 


DO O07 L=1,50 
DOWN= (CON*TAW*EXP (ROW*X (L) )) + (CCON*SI/EXP (+ROW*X (L) )) 
Q +((PI*MEW+ (BETA/SIGMA1)) /EXP (+SIGMA*X (L) )) 
UP=(CON*SI*EXP (ROW*X {L) )) + (CCON*TAW/EXP (+ROW*X (L) ) ) 
Q +#((PI¥MMEW-(BETA/SIGMA1) ) /EXP (+SIGMA*X (L))) 
ABS (L,2Z) =( (PVK1* (DOWN+UP) ) + ((IN*(1-RC) ) /EXP 
Q (SIGMA*X (L))) 
Q *PVK) *DELX 
IF(L.£Q.1) BCT=(RC*IN) + (1- (RC¥CTI) ) * ( (CCON*TAW) 
Q +(CON*SI) + ((PIXHMEW- (BETA/SIGMA1) ))) 
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IF (L.£0.50) BCB=(1-(RC1*CTI) ) * ( (CON*TAW*EXP (ROW 
Q ¥*TH(LL))) + (CCON*SI/EXP (ROW*TH (LL))) + 
Q ((PI¥MEW+ (BETA/SIGMA1) ) /EXP ( 
Q SIGMA*TH (LL))))+( (IN* (1-RC)) /EXP (SIGMA*TH (LL) )) 
0 *(1-RC1) 
CONTINUE 

GO TO 08 

CALL ISHEET (ROW, SIGMA, PI,CI,DELX, PVR, PVK 
Q ,RC,SIGMA1, BETA, BCT) 

DO 09 J=1,50 

ABS (J,Z) =CABS (J) 

CONTINUE 

GO TO 98 

TH1=TH (LL) 

CALL ABSORB (DELZ EVR, RC, SIGMA, THI, BCT, BCS) 

po 14 7JJ=1,590 

ABS (JJ, Z) =CABS (JJ) 

CONTINUE 

CALL RAD BAL (BCT,BCB,Z,DELWL) 

CONTINUE 


pO110 KK=1,27 
TOTIN=0.5* (TOT (KK) +TOT (KK+1)) #+TOTIN 
CONTINUE 


TOTAL ATTENUATION 
DO oi, L=1950 
O0e12az=Ie2s 
BABS=(DELWL*9.5) * (ABS (L, Z) +ABS (L, Z+1)) 
IF (Z.EO.1) TABS (L)=0.0 
TABS (1) =BABS+TABS (L) 
CONTINUE 
TABS (L) =TABS(L) /TOTIN 
ATTEN (LL) =TABS (L) ATTEN (LL) 
CONTINUE 
BCT=BCT/TOTIN 
BCB=BCB/TOTIN 
TRANS (LL) =BCB 
REFLEC (LL) =BCT 


Hog1s9rre1950 

EXC (II) =TABS (II) / (TH (LL) /50) 
CONTINUE 
TEXC=ATTEN (LL) / (TH (LL) ) 
SUM=ATTEN (LL) +BCB+BCT 
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17 


WRITE (6,200) 
WRITE (6,1100) T,TH (LL) 
WRITE (6, Fao) PYR 
WRITE (6,599) 


WRITE (6,409) (X(L) , TABS (L) , EXC (L) , L=1,50) 


WRITE (6,300) ATTEN (LL) 
WRITE (6,300) TEXC 
WRITE (6,310) SUM 
WRITE(6,1000) BCT, BCB 
CONTINUE 

RETURN 

END 


ENERGY---CALCULATES INCIDENT RADIATION ASSUMING BLACK 


BODY DISTRIBUTION, 


SUBROUTINE ENERGY 


REAL IN 
COMMON WL,IN,T,X (50) ,CABS (50) 


ee 4387) /(#L*T) ) 


ET= (EXP (ET) -1.0) 


ae ( (3. 7498-05) /( (HL¥*5) *VET) ) /3. 14159 
RETURN 


END 


ABSORB-~--DETERMINE ABSORBED, REFLECTED, AND TRANSMITTED 


RADIATION OCCURING WHEN SCATTERING COEFFICIENT 
IS ZERO OR INSIGNIPICANT COMPARED TO THE 
ABSORBTION COEFFICIENT. 


SUBROUTINE ABSORB(DELX, PVK, RC, SIGMA, TH1, BCT, BCB) 
REAL IN 

COMMON WL,IN,T,X(50) ,CABS (50) 

DONO 1 edad950 

IF ((X (J) *PVK).GT.50) GO TO 03 

ABS 1= ((1-RC) ¥ IN) EXP (-PVK* (X (J) - (DELX/2)) ) 
ABS 2= ( (1-RC) *IN) EXP (-PVK* (X (J) + (DELX/2) ) ) 
CABS (J) =ABS1-ABS2 

IP (J.EQ.1) BCT=RC*IN 

IF (3.20.59) BCB=(IN*(1-RC)) /EXP (SIGMA*TH1) 
GO TO 01 

CABS (J) =0.9 

IF (J.EQ.1) CABS (1) =(1-RC) *IN 

IP(J.EQ.1) BCT=RC¥IN 

IF (J.EQ.50) BCB=0.0 

CONTINUE 

RETURN 

END 
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BAL---SUMMATION OF ALL RADIATION ABSORBED,REFLECTED 
,AND TRANSMITTED, 


SUBROUTINE RAD BAL (BCT,BCB,Z,DELWL) 
INTEGER Z 

DIMENSION 1(28),B(28) 

T (Z)=BCT 

B(Z) =BCB 

IR (Z7LTA2a8y GO TO 01 

CB=9.9 

CT=0.0 

DORG29I=1, 27 

TT= (DELWL*9.5) * (T(J) +7 (J+1)) 
BB= (DELWL*0.5) * (B(J) +B (J+1)) 
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EXP (SIGMA*X (J)) 
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IF(J.EQ.1) BCT=(RC*#IN) + (1-(RC*CTI) ) *(((R/B) 

*TAW) + ((PI¥MMEW- (BETA/SIGMA1)))) 

60 10 104 

TE((SICHA¥ XJ) ) GE. 70)) GO Tons. 

CABS (J) =(( (IN* (1-RC) ) /EXP (SIGMA*X (J))) 
*PYK) *DELX 

IF(J.EQ.1) BCT=(RC*IN) 

ef) ae) 8 

CABS (J) =0.9 

IF(J.EQ.1) BCT=(RC*IN) 
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RETURN 

END 
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Type of 
Equipment Used 


Blackbody Source 


Chart Recorder 


Laser 


Radiation Probes 
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Description 


In order to calibrate the thermal 
radiometer, a known intensity of 
radiation was produced by a Black- 
body Radiation Source Model WS143 
made by Electro Optical Industries, 
Inc. The selected source tempera- 
tures were maintained by a Tempera- 
ture Controller Model 215B pro- 
duced by the same company. 

A Hewlett Packard strip chart 
recorder Model 7100B was used on 
the 0.5 millivolt and 1.0 milli- 
volt scales to record thermo- 
couple outputs. 

The specific laser used was a 
Spectra PhysiicsmModelelsOB. eeihis 
laser was of helium-neon-gas 

phase type producing continuous 
emmission at a wavelength of 

ORGS 2a: 

To measure radiation intensities 
two types of radiometers were 
used. They are a Cintra Quantum 
Radiometer Model 101 and a Cintra 


Thermal Radiometer Model 202. 
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Gr. 


Type of 
Equipment Used 


Radiation Sources 


Thermocouple 


Calibration 
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Description 


To detect the radiation intensi- 
ties the quantum radiometer used 
a Cintra Probe #1114 while the 
thermal radiometer used a Cintra 
Probe #2000. 

The tungsten lamp (3000°K) con- 
Sisted of a Sylvania Tungsten 
Halogen Lamp type DVG rated at 
650 watts surrounded by a para- 
Doli Cane tlectone. = Omsimulate ache 
solar spectrum, a Metro-Lite 
Model ME-4 carbon arc lamp was 
employed. 

To calibrate the thermocouples, 

a platinum resistance thermometer 
cat. #8163-QB produced by Leeds 
and Northrup was employed. Out- 
puts from the copper constantan 
thermocouples and resistance 
thermometer, that were obtained 
by placing them in a temperature 
controlled bath, are compared using 
Potentiometer #7555 Type K-5 and 
a Guarded Nanovolt Detector 
#9383-1 also produced by Leeds 


and Northrup. 


oct at hl, AEE Asa T RMD S 
srditd b bedu vodonotben Vemamdy 
| ee 
enon (#°@D02) gael aateginsd ent 
asdzonuTt sinevtye & TH wstete 

be beset 990 eqys qmad eapetet 
~st6q 6 Ud -bebewor ie 2aoeWweeD 
a3. esn lente oT otae fren Ot ies 
attJ-cuseh o -ntycoresye doing 
Pow-qne! are hodveo F-24 We00N 
abeybt gif 

, eet qyuoonmrend oxft sdanvita of 
sorgmont sls. s20R22he9% nombre hy 5 
bees Ys bemwthowg SP-E ara sted 
-tu0 bevel qe crew queers tae: 
memngtencs ASggns ent AON etiG 


ca 
“aa ai 
veh iy 


